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ABSTRACT 

The basic idea of baryogenesis is lectured to introduce non-experts to this subject. 
Some recent topics, necessarily subjective in view of short time limitation, are also 
presented to show how the initial condition for baryogenesis is realized in the new 
framework of inflation. 
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1 Introduction 



Let us first recall the modern view of cosmology; big bang theory. According to 
general relativity the spacetime evolution is determined via the Einstein equation 
by the matter content of universe, which differs from epoch to epoch depending 
on what kind of energy dominates the energy density of the universe at that time. 
There are three important epochs characterized by different pressure-energy relation; 
(1) vacuum energy dominance with p = — p , (2) massless particle dominance with 
p = ~p , (3) nonrelativistic particle dominance with p <C p . With homogeneity 
and isotropy of 3-space the Einstein equation is much simplified with the Robertson- 
Walker metric: ds 2 = dt 2 — a 2 {t) x 2 . The behavior of the scale factor a(t) then 
follows; 

(1) a(xe H \ H=^MGV/i, (1) 

(2) a oc t 1/2 , (2) 

(3) a oc t 2/3 , (3) 

in their respective energy ranges. The first stage is the inflationary epoch where a 
constant vacuum energy V gives the exponential growth of the scale factor, which 
is believed to solve the horizon and the flatness problems in the old big bang model 

It would be of some interest to view the Hubble parameter defined by H = a/a 
as a function of a or a typical energy scale T oc 1 ja of dominant massless particles in 
the universe, instead of a(t) as a function of time. In this picture the first inflationary 
stage is longest, lasting for the scale change > 10 30 . Subsequent radiation and matter 
dominant epochs have energy changes of order, 10 26 or 10 3 . The standard model of 
microphysics only probes the temperature range of 10 2 GeV — 1CT 13 GeV , of much 
smaller variation than in the inflationary epoch. Of great importance in subsequent 
cosmic evolution and of intensive current interest is the transient stage from inflation 
to radiation dominance; the epoch of reheating after inflation, which I shall touch 
upon later in this lecture. 

For a long time the success of the big bang model rested with three cornerstones; 

(1) Hubble expansion, 

(2) Planck distribution of relic photons, 

(3) light element abundance such as 4 He, D, 7 Li. 
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In the last decade there have been substantial improvements of this classic achieve- 
ment both in observation and in theory. Perhaps the most significant is COBE mea- 
surement of the spectrum shape and large scale fluctuation of the relic microwave 
temperature: absolute motion of earth relative to the universal cosmic expansion has 
been observed with 5T/T rs 1CT 3 , and after subtraction of this dipole component 
there exists a large scale fluctuation || of order 10~ 6 , whose precise feature yet to be 
clarified is very important in the theory of structure formation. Nucleosynthesis || 
has become a precise science developed to a fine detail, so much refined that one can 
even discuss dependence on the number of neutrino species and the neutron lifetime. 
But presumably the most important in nucleosynthesis is a consistent and reliable 
determination of the key cosmological parameter, the baryon to photon ratio Ub/tIj, 
of order 10 -10 . 

In order to probe microphysical processes that occur at each instant of cosmolog- 
ical time, it is necessary to compare the two time scales; physical time scale during 
which a relevant process takes place, and the Hubble time scale of order 1/H. If 
the time scale of a physical process is much shorter than the Hubble time, then this 
physical process will occur frequently, and in the limit of large rate equilibrium is 
reached. 

In a weak coupling theory such as supersymmetric grand unified theories beyond 
the TeV scale, the proper framework to discuss the interaction rate is the Boltzmann 
equation for one particle distribution function taking into account the cosmological 
expansion. When the Boltzmann equation is integrated with respect to the momen- 
tum phase space, an equation for the number density follows; 



where r,[/] is the invariant interaction kernel that depends on all distribution func- 
tions fj participating in the process. As an example let me give the form of the 
kernel when only two-body reactions are involved, 



where ^ij^ki is the invariant rate for the process ij — > kl essentially given by the 
differential cross section times the initial relative velocity, and ± corresponds to 




(4) 



Uf] = \ll I dVjdVkdVi [ (1 ± /;)(! ± fj)fkfi in 



j ,k ,i 



-(l±f k )(l±f l )f i f jlij ^ kl ], 



(5) 
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either stimulated emission for bosons or Pauli blocking for fermions. dV 



dC A p 



is 



£(2tt) 3 



the invariant phase space element. 

In general it is difficult to solve non-linear integro-differential equations of the 
Boltzmann type, but it is often possible to estimate reaction rate in the right hand 
side by some means. In certain situations the rate may roughly vary linearly with 
the number density itself nf, 



with H — a/a. For instance, if the relevant particle freely decays having no inter- 
action with medium particles and no process of particle production operates, is 
simply the free decay rate. 

A particular reaction is of cosmological relevance only if the reaction rate given by 
the right hand side of the Boltzmann equation is much larger than the Hubble rate 
3Hrii. What typically happens is that the process continues at early times and after 
some late time it is frozen. In cosmology it is thus important to estimate decoupling 
epoch of a particular reaction. The crucial quantity in determining relevance of the 
reaction is the ratio of the average rate to the Hubble rate H. If this ratio is very 
large, the reaction occurs frequently. 

As an important practical example let us consider the four Fermi type interaction 
governed by a coupling a/M 2 . This could be the weak interaction involving neutrinos 
such as vv — > e + e~ or baryon number violating reaction, qq — > ql , if one uses as 
the relevant mass M the electroweak or the GUT mass scale. Without any detailed 
discussion I would give you a suggestive form of the average rate, 



Here T is a typical energy of participating particles and of order the temperature 
if all these particles are in thermal and chemical equilibrium, and the last factor 
T 3 represents the number density of colliding particles. The factor in front of this 
number density is the averaged reaction rate for the elementary process, which grows 
with energy T as a 2 T 2 /M 4 for T <C M, but is finally saturated and decreases as 
a 2 /T 2 in the high energy limit. On the other hand, the Hubble rate is of order 
v / iVT 2 /m pl with m pl = 1/y/G ~ 10 19 GeV the Planck mass. N is the number of 
massless particle species in equilibrium. It is then easy to see that the ratio 




(6) 



(T 2 + M 2 ) 




(7) 
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has a maximum at around T = M and the maximum ratio is of order 

a 2 m pl /(VNM) . (8) 



The condition of the large rate for equilibrium is thus M <C a 2 m p \/ 'yN . Since 
the right hand side contains the Planck mass, this condition is easily obeyed for the 
weak processes. This is an important observation to understand why the element 
abundance cooked up in the early universe is insensitive to the initial condition: one 
can justify that the initial neutron-proton ratio is given by the value in thermal and 
chemical equilibrium. 

As another example let us consider how much of baryons is left over if the uni- 
verse is symmetrical with respect to baryons and anti-baryons. This is a problem 
relevant to baryogenesis, the main theme of this lecture. As the final element car- 
rying the negative baryon number, anti-nucleons gradually pair-annihilate via many 
pion production; NN — > many tc' s , but after the decoupling of this annihilation 
process the leftover abundance n^jn^ = n^/n 7 is frozen. The annihilation process 
has a cross section roughly of order (av) wl/ m 2 at low energies, independent of the 
energy participating in the process. On the other hand, the thermal number den- 
sity is suppressed by the Boltzmann factor when nucleons and anti-nucleons become 
nonrelativistic as actually is the case; 

n N = n w ~4(^r/ 2 e- m "/ T . (9) 

From the decoupling condition, (av) • un = H , one then finds for the decoupling 
temperature T d and the leftover abundance that 

T d ^ [ln m N m pl (av) } _, ^ 1 (1Q) 



m N m^/N 45 ' 

1(T 18 . (11) 



n N HlvA - im _ 1s 



n 7 mjqm p \(av) 

The final abundance un/uj is too small compared with the observed value ~ 10~ 10 , 
and this argument essentially rules out the symmetric cosmology. 

In the asymmetric universe with respect to baryons and anti-baryons the present 
number ratio %/n 7 is the measure of the imbalance between baryons and anti- 
baryons prior to the annihilation process: 

M , (12) 

V ^7 / present V ) before annihilation 
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since in thermal and chemical equilibrium baryons, anti-baryons, and photons are 
all in roughly equal abundance, B « B m 7 , ignoring the mass threshold and the 
statistical factor. Our task is thus to explain this quantity from microphysics. 



2 Condition for baryogenesis 

Although necessary ingredients for the baryogenesis were written a long time ago[Q, 
intricacy of this condition has been spelled out much later [fj,0,[0,||. (For reviews, 
see [H,[10|.) Their fascinating features still invite many interesting scenarios. The 



three obvious conditions are 

(1) baryon number non-conservation, 

(2) violation of dicrete C and CP symmetry, 

(3) departure from thermal equilibrium. 

The baryon number is not a sacred symmetry in modern gauge theories, and 
indeed violated explicitly by heavy X boson mediated processes in grand unified 
theories |TT|. Moreover, even in the standard SU(3) x SU(2) x U(l) gauge theory 
the baryon number is violated by instanton effects at zero temperature ||12| , which is 
however unobservably suppressed. It has however been recognized that sphaleron- 
mediated processes [13| may enhance baryon non-conservation in the standard theory 



at high enough temperatures, higher than TeV scale, although B — L is exactly 
conserved here. The sphaleron |TJ]] is a finite energy field configuration which bridges 



between different nonperturbative vacuum configurations ordinarily suppressed by 
tunneling probability at zero temperature. The effect may be enhanced at a high 
temperature T by the factor e~ F ' T , with F the free energy of this unstable sphaleron. 
The sphaleron effect can readily wash out the baryon asymmetry generated prior to 
this epoch if initially B — L = 0. It should however be noted that it is not a trivial 
matter to create from a symmetric state the asymmetry at the electroweak scale. 
This is related to that the last two conditions are much more subtle, as I shall 
explain shortly. 

Is the standard model capable of explaining the baryon asymmetry? First of all, 
the first two conditions are met as a matter of principle even in the standard model, 
although strength of B violation is to be detailed. The last condition of departure 



from equilibrium is in general difficult to meet. Most scenarios [15] use the first order 



phase transition of the electroweak gauge symmetry breaking as a means of setting 
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up non-equilibrium environment. This requires a Higgs boson mass < 45 GeV in the 
standard model of one Higgs doublet, which already seems to be ruled out by LEP 
experiments. But this bound is based on a one-loop computation, so there might be 
some loophole in the argument not contemplated so far, and it may be worthwhile 
to explore the possibility of baryogenesis, ignoring the Higgs mass bound. 

Farrar and Shaposhnikov [|l!J recently proposed an interesting scenario that em- 
ploys quark scattering off the electoweak bubble created at the first order phase 
transition. As an idealization one may assume that baryon number is strongly vio- 
lated in the unbroken phase, while it is conserved in the broken phase. Quarks are 
scattered off within some length near the bubble wall and there may be many scat- 
tering amplitudes that interfere, giving non-trivial CV violating effect as required. 
Unfortunately the original calculation neglected the important effect of coherence in 
the cosmic plasma [|17]. Coherence is crucial to generate the baryon asymmetry by 
this mechanism and is only maintained by the coherence length which is severely 
limited in the cosmic plasma. The resultant asymmetry is too small to yield the 
one needed by observation, typically < 1CT 20 ]T7|]. Thus at the moment there is no 
viable model of baryogenesis that explains the observed ratio using only the source 
of CV violation identical to the CKM matrix. 

According to our present wisdom the standard model must be extended in order 
to explain the baryon asymmetry. There are two attitudes in extending the estab- 
lished physics with this respect: one is to use the new physics of unification such as 
grand unified theories or its supersymmetric extention that have the natural source 
of all ingredients for baryogenesis, and the other |L5| is to minimize the extention 
utilizing the strong electroweak baryon violation maximally. The way I would like to 
characterize these two directions is that the first direction needs a big jump beyond 
the standard theory, while its physical mechanism of baryogenesis is straightforward 
and readily understandable. On the other hand, the second extention needs a small 
step beyond what we already know, such as a slight modification of the Higgs system, 
but physics involved is fairly complicated, and there seems no model everyone would 
like to consider seriously. With this situation in mind I would like to mainly discuss 
essential features of baryogensis based on GUT. I shall also discuss obstacle against 
the GUT scenario and how to evade it. 

Although not discussed below, there is another interesting scenario of baryo- 



genesis that is becoming popular recently; the Affleck-Dine mechanism [13. This 
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mechanism uses a feature of supersymmetric models; existence of many flat field 
directions in the scalar field potential. These directions may include field condensate 
that carrys the baryon number or the lepton number, both of which can be used 
when combined with the strong electroweak baryon nonconservation at finite tem- 
peratures. The resulting baryon to photon ratio tends to be large compared to the 
observed one, and one usually has to consider some process of dilution at the same 
time. A version of this mechanism is discussed in the lecture by L. Randall at this 
school. 



3 GUT baryogenesis 



Baryon nonconserving processes in grand unified theories are mediated by gauge or 
Higgs bosons generically called X bosons. The simplest and workable scenario of 
GUT baryogenesis employs X boson decay, which has two types of decay modes, 
the two quark decay mode qq, and the leptoquark mode ql. Coexistence of the two 
modes with different baryon numbers is a manifestation of baryon nonconservation. 

Two simple examples of grand unified theories are SU(5) and SO(10) models 
fill . These models have unified multiplet structure of one family of quarks and 
leptons of the form, 



£77(5); 5 + 10, 

( d R \ ( Wb 
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V v ) 

50(10); 16, 
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u R u g u b N 
d R d G d B e 



(13) 



(14) 



For instance, in the £77(5) model X gauge bosons induce transitions from a quark 
to a lepton in the representation 5, while in 10 they cause transitions both from a 
quark to a lepton and from an anti-quark to a quark. 

Recent precise measurements of coupling constants at LEP suggest |L9| that 
supersymmetric extention of the SU (5) model gives a consistent picture of coupling 
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unification using the renormalization group. The 5*0(10) model, on the other hand, 
is interesting if a finite neutrino mass hinted by the solar neutrino experiments is 
real. Despite of the lack of evidence for proton decay grand unified theories thus 
deserve serious consideration. 

In discussing the GUT baryogenesis, I shall dismiss complication due to many 
decay channels and concentrate on the two decay modes for simplicity. Let us denote 
decay rates of the two modes by 7 9 , 7/ for X and 7g , 72 for its anti-particle X. 
When a pair of X and X decays, a finite baryon number may be created with a rate, 

2 12 1 

AB = 3 7<? ~ 3 7i ~~ 3^ q + 3 7z = lq ~ 7<? ' ^ 

This quantity was simplified by using the requirement of CPT theorem, which states 
the equality of particle and anti-particle total decay rate; 

Iq+ll =lq + ll- (16) 

Clearly CP violation is called for 7 g 7^ 

Before discussing perturbative calculation of this asymmetry, I shall explain the 
non-equilibrium condition for the GUT baryogenesis. Departure from equilibrium 
demands in this case that the inverse decay process, qq — > X , ql — > X etc. 
is essentially frozen by threshold effect. This is realized if at the decay time the 
temperature is too low to create heavy X's: T < mx when the Hubble H(T) = 7^ . 
This leads to O [y/m v \ jx } < m x ■ Since the total decay rate 7x ~ am x , it gives a 
constraint on the X mass, 

m x > 0[am p \] . (17) 

A more precise estimate yields the X mass bound of order 10 15 — 10 16 GeV, close to 
the GUT scale. This constraint may be obeyed without much difficulty by a Higgs 
X boson, although with some difficulty by the gauge X boson. 

Calculation of the magnitude of the baryon asymmetry due to the X boson 
decay involves interesting interference effect. Suppose that one computes the baryon 
production rate AB in perturbation theory; 

AB = ]T \ gi f l +g 2 f 2 + ...f-\glf 1 +g* 2 f 2 + ...f 

phase space 

= -4^(g l9 *)- £ + • • • ■ (18) 

phase space 
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Each decay amplitude gifa corresponds to a Feynman diagram and gi lumps together 
all coupling factors leaving the dynamical part of the amplitude to fi. The rate for X 
thus has the complex-conjugated coupling g*, as required by CVT or hermiticity of 
the Hamiltonian. The final expression of AB clearly indicates need for a dynamical 
phase ^{fifZ) 7^ besides the CV violation, Q^gig^) 7^ . The dynamical phase may 
arise as a rescattering phase in ordinary two-body processes, and one may regard 
the phase above as a generalization of the rescattering phase. The constraint for a 
non-trivial dynamical phase is that intermediate states must have a threshold below 
the parent particle mass nix- This observation also hints an efficient way to compute 
the asymmetry as a discontinuity according to the Landau-Cutkovsky rule: one can 
put intermediate particles on the mass shell which makes actual computation of the 
asymmetry much easier. 

The baryon to photon ratio, or more precisely the baryon to entropy ratio %/s 
may be derived from the baryon production rate AB as follows. Define first the 
baryon number in the comoving volume, Mb = a 3 ^s, together with a similar X 
density, Mx = a 3 n X - These vary as 

Mb = AB ■ Mx , Mx = - 7tot ■ Mx , (19) 

with 7 tot = 7 g + 7z the total X decay rate. This is integrated with the condition that 
initially no baryon number exists and finally no X boson exists: 

A R 

(M B )f = e(M x ) l , e= . (20) 

7tot 

Since this reaction and subsequent evolution proceeds almost adiabatically, the en- 
tropy per comoving volume is approximately conserved; a 3 s = constant. This 
relation can be used to eliminate the volume factor a 3 in favor of the entropy density 

(^)/ = e(^) ( . (21) 
s s 

There are many interesting details of how the magnitude of the baryon asymmetry 
may be correlated with some other physical quantity in a particular model. For 
instance, the 577(5) model with the minimal Higgs structure yields too small a value 
of Ub/uj. On the other hand, the baryon asymmetry is correlated with the neutrino 
mass in 50(10) models [pCfl : t ne smaller the neutrino mass is, the smaller the baryon 
asymmetry is. This can be utilized to constrain a finite, but a small neutrino mass. 
I shall however omit discussion of these subjects in this introductory lecture. 
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It is now appropriate to discuss difficulty associated with the GUT baryogen- 
esis and how to overcome it. First, as I already mentioned above, there may be 
baryon number annihilation via electroweak processes at finite temperatures. One 
can always excuse for the GUT scenario that the electroweak baryon nonconserva- 
tion keeps intact B — L, so if B — L is created at the GUT epoch, later evolution 
only redistributes B and L, but never annihilates the baryon number; 

-Bfinal — C - (B — L) initial , (22) 

with Lfinai = (1 — c) ■ (B — L)i n itiai • Here c is a calculable number of order unity. 

One has to be careful however to avoid a large lepton violation at intermediate 
temperature scales that may potentially dissipate away the baryon number with 
AB = AL = , when combined with electroweak B nonconservation. Allowed 
range of this interaction, accordingly constraint on the neutrino mass via the seesaw 
type of neutrino mass generation may be estimated by introducing a generic type of 
AL = 2 interaction; 

£al^o = Ihhw + (h.c.) , (23) 

with l L and (p the lepton and the Higgs doublet, respectively and v ~ 250 GeV. 
Requirement of harmless lepton violation, Tal^o < H , leads to the upper bound 
on the neutrino mass f2"Tf ; 

m u < 4 x KT 3 eV (T B _ L /10 16 GeV)~ 1/2 , (24) 

with Tb-l the scale of B — L generation at higher temperatures. This constraint 
is always applied to the lightest neutrino species, but when the neutrino mixing is 
large, it is also applied to the heaviest neutrino species. 

The second problem with the GUT baryogenesis is more serious; a possible 
overproduction of gravitino and associated low reheating temperature after infla- 
tion. The gravitino is the superpartner of graviton in supergravity theories. It is 
a spin 3/2 particle and couples with ordinary matter field with the gravitational 
strength. Moreover, one usually associates supersymmetry breaking scale with the 
electroweak scale in order to ease the hierarchy problem, thus the mass of the grav- 
itino m 3 / 2 = 0[TeV]. This has a consequence potentially very serious, because the 
decay rate of the gravitino is given by 

r^lJ^secrM^, (25) 
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the lifetime being close to the epoch of nucleosynthesis. If the gravitino abundance 
is larger than of order 10~ 10 relative to the entropy density, the successful nucleosyn- 
thesis would be destroyed. 

For a short while inflation was considered to save this potential disaster by di- 



luting away the gravitino abundance. But subsequently it has been recognized [22 



that regeneration of gravitinos after inflation severely constrains the maximally al- 
lowed reheating temperature. The point of this argument is that gravitino pairs may 
be produced from ordinary particles, whose cross section is calculable and of order 
l/jripi, giving the abundance of order 

^ = O[10- 2 ] — , (26) 
s m p i 

in terms of the reheat temperature T after inflation. For successful nucleosynthesis 
one gets the constraint on the reheat temperature 

T<O[10 10 -10 n ]GeV (27) 

from < O[10~ 10 — 10~ n ] . This temperature is too low to create the baryon 
asymmetry by the X boson decay. It thus appeared that supergravity models arising 
as the field theory limit of presumably the ultimate theory of superstring give too 
low a reheat temperature incompatible with the GUT baryogenesis. 

I shall discuss in the following sections how a correct theory of reheating after 
inflation may provide a high temperature phase suitable to the baryon generation 
without overproduction of gravitinos. 



4 Thermal history after inflation 

How the hot big bang is started after inflation is a fascinating subject that can be 
discussed independently of baryogenesis. But since the subject directly addresses the 
origin of entropy in our universe at earliest times of evolution, one may very natu- 
rally entertain the possibility that the two basic quantities in cosmology, the baryon 
number and the entropy are both created roughly at the same time. Moreover, as- 
sociation of the origin of the cosmic entropy with inflation sets an ideal theoretical 
framework, because inflation dilutes away everything in our observable part of the 
universe: one must explain the origin of entropy starting from empty space, except 
the coherent inflaton oscillation around the minimum of inflaton potential. 
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For quite some time the theory of particle production due to coherent and homo- 
geneous field oscillation was based on a naive picture of lowest order perturbation 
p3| . The oscillating inflaton field, denoted here by £(£), is regarded in this naive 



picture as an aggregate of condensed bosons at rest, and these bosons are assumed 
to decay stochastically according to the rate given by the Born approximation. Thus 
equating the decay rate with the Hubble rate 7^ = H , together with assumption of 
the instantaneous reheating, leads to 

T B ~ 0.1 ^m-[ ~ 10 14 GeV ■ 9 ■ ^ ^ . (28) 

The Born rate for massless particle decay, £ — > ipip , 

■<^w^ (29) 

was used in this estimate. 

This picture is valid when the amplitude of oscillation and the coupling g to 
matter field is small enough, but is grossly wrong if this condition is violated, for 
instance when the amplitude of oscillation is large. As will be explained shortly, there 
exist an infinitely many bands of instability that may contribute to the inflaton decay, 
only one of which, the lowest band, when restricted to the small amplitude limit, is 
identified with the Born decay rate. This phenomenon of instability is known as the 
parametric resonance under periodic potential plj , |25[ . 



A systematic method to understand particle production and associated inflaton 
decay can be formulated PB|,||27[| in the Schrodinger picture of quantum field the- 



ory . (For other approaches and other aspects of this problem, see ref ||28||,||29||,|]30 



Quantum bose fields that couple to the inflaton field are treated as a quantum 



operator (in the Heisenberg picture), but the inflaton field is regarded as classical 
in this approach, although back reaction against the inflaton oscillation is also con- 
sidered. In the Schrodinger picture the state vector describing behavior of quantum 
field coupled to the inflaton is given by a direct product of state vectors \ip(t))^ of 
independent spatial Fourier modes. During the time interval At that obeys 

l/m € <At<l/#, (30) 

one may assume exactly periodic inflaton oscillation with a periodically varying 
frequency, , containing the oscillating function £(£). For longer time scales, 

At ^> 1/H , the amplitude of oscillation is taken to adiabatically change. 
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A salient feature of this approach is that one may solve the quantum state with 
a Gaussian ansatz; 

(g*W%=-^=exp[I^gJ], (31) 

y/Uk 2ll k 

where u k (t) is shown to obey the classical oscillator equation; 

[ ^ + wf(t) ]u k = 0. (32) 

The initial condition for this classical equation is determined by the choice of an 
initial quantum state. The simplest, and a reasonable choice is the ground state 
with respect to some reference frequency, most naturally the frequency at the onset 
of inflaton oscillation 00^(0). This is because after inflation the state is essentially 
devoid of matter. Under this circumstance the initial condition is 

u fc (0) = (^)- 1 / 2 , ^(0) = z^ , (33) 

71 U k 

with uq the reference frequency, and one may further simplify the wave function in 
terms of |"U fc (t)| alone p6[| . 

The most important consequence PE| ,||27|| of this formalism is that it gives a 
rationale to introduce a coarse grained density matrix p^ which has a classical 
probability distribution. The coarse graining here is defined by a short time average 
over the time scale of order a few oscillation periods. This seems a reasonable way to 
extract global behavior of the quantum system ignoring fine details of the quantum 
state. This makes it possible to replace the quantum density matrix by the time- 
averaged diagonal part (n\p^\n) of the density matrix in a convenient base such 
as the Fock base of frequency u> : 

(n\pM\n) = (n\m)(m\n), (34) 

where \n) is the n— th level of field oscillators and the overline represents the short 
time average. After the coarse graining a finite entropy may be assigned; 

- tip {D) \np {D) > 0. (35) 

From this reduced density matrix one computes various physical quantities. For 
instance, the produced particle number in each mode is given by 

1 - 1 J! , 1,2 „2X _ 1 \ _ tv /J fl ).t. 



(N k ) = {-^pi + ^qi)--)=tT{ P ra\a k ) 



4( uj\u k 


2 -7T) 


2 + 


( d\u k \ 2 y 


16tiuj 




2 



(36) 
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with q k , p k oscillator coordinates. Suitable time average is understood here. As 
t — > oo, 

(N k ) -> e Afem «* x (polynomial in i) (37) 

for the momentum k in the instability band. 

Two types of matter coupling have been considered; 

V^/CV + ^c*?™^ 2 , (38) 
with c a constant of order unity and the mass of the inflaton. As a model 



of inflation we consider the simplest chaotic type inflation |32j with parameters; 

~ 10 13 GeV consistent with COBE anisotropy, and the initial amplitude of 
order the Planck scale, 

£o~yJ™ P i- ( 39 ) 

The dimensionless coupling g is taken arbitrary at the moment. 
The classical mode equation is then 

d 2 u 

— + [h-29 cos(4z) - 4cVe sm{2z) \u = , (40) 

h ^ 4 ^L + 2 e, * = , = 3*. (4i) 

where is the amplitude of inflaton oscillation; 

£(t) = -£ sin(m € t). (42) 

The criterion of large or small amplitude is thus given by the magnitude of 9. In the 
very small (9 <C 1) or in the very large (9 ^> 1) amplitude limit the classical mode 
equation ( fiO|) reduces to the Mathieu equation |33| with a single oscillating term. In 



both these limits we developed analytic formulas suitable for detailed analysis of the 
reheating problem and related problems, too [EfJ,[33 . 



The structure of instability bands is as follows. Each band is labeled by an integer 
n = 1 , 2 , 3 ■ • •, and goes to h — > n 2 in the small amplitude limit, £o ~~ > 0. The band 
width in the small amplitude region is 

6> n/2 c 

AK = 2^[(n-l)!] 2 ' (43) 



at a fixed 9, with c n some function of c |[26|,|[35||. For instance, c n = 2™ with the 
Yukawa coupling | gm^ip 2 alone. The instability bands are thus indeed very narrow 
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in the small amplitude region. But as 9 increases, the bands become broad, and for 
h < 29 most of the (9 , h) parameter space is covered by instability bands except 
bounding narrow stability bands. For parameters (h , 9) or corresponding (k , £o) 
within an instability band the classical solution exhibits exponential growth; 

u -> e Am «* x P(t) , (44) 

with P(t) a periodic function. This has the important consequence of exponential 
particle production rate unless the back reaction stops it. 

In the small amplitude limit exact results for the growth rate A and mode sum 
are available. After the coarse graining the initial state decays according to e~ rvt 
where 

r = £r„, r n = ^ £ x s . (45) 

n k e n-th band 



T n is the decay rate per unit volume of the n— th band, which is computed 

as 



4 

r - 5 



1 - 4^ (Ah n ) 2 . (46) 



64vr \ 

It was also shown recently that the small amplitude result can be understood 
by familiar perturbation theory. The point is that the mode-summed rate grows 
with the amplitude as oc £q™ oc {n^) n where = | ^ is the number density 
of condensed inflatons. The decay rate precisely coincides with the zero-momentum 
limit of n to 2 body process n £ — > (ftp , worked out using the ordinary Feynman 
rule. In particular, the decay rate of the first band with c = 1 is given by 

,2^,2 c2 



1 64vr 



\ 



K . (47) 



This exactly coincides with the one-particle decay rate 75 for £ — * when divided 
by the inflaton number density \ m^^Q. 

It is important for many interesting applications to work out analytic formulas 
in the large 9 region. This seems a fomidable task in view of that non-perturbative 
analysis is involved. A remarkable result |M| is that in the functional Schrodinger 
picture one can solve the fundamental quantity Uk{t), which becomes rigorous in the 
region, 9 ^> 1 with \h — 29\ 9 . Furthermore this is precisely the parameter region 
of prime importance when one considers cosmological evolution, as will be discussed 
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shortly. In the rest of this lecture we shall neglect the mass m of quantum bose field 
considering only those of m <C rn^. 

Cosmological evolution changes both the momentum and the inflaton amplitude 
according to 

k (x l/a , f oc l/a 3/2 . (48) 

Taking in the large 9 region the leading variation alone gives the approximate rule; 
Ah ~ 2 A9 . Thus parameters of the dominant contribution move parallel to the 
h = 29 line with cosmological evolution. Moreover, the largest particle production 
with the largest rate A occurs in the deepest region within instability bands of the 
Mathieu equation (the original two terms of oscillations reduce to one term for 
9 ^> 1). This implies that the most dominant region is along h = 29 within some 
width 5 = 5{h-29). 



The analytic result derived in ref is summarized as 
1 



A = - ln(Vx + Va 73 !) , (49) 



7T 



x = (l + e -^/( 4 ^)) cos 2 V, (50) 



TV 2 



^ = !LV0 + ^ln(^V)+-31n[-^ '( ]. (51) 

4 4V9 V V2 ' 2 [ r(| + ^/(8v^)) 



The instability region is characterized by x > 1, while the stability region by x < 1. 
These two regions alternate roughly with equal band width of A9 ~ around 
along h = 29. For 5 or \J 9 m ^0i it approximately follows that 

1 7T 2 

A = — ln( v 7 ^ | cos ^ I + V I cos(2^) | ) , i/j = — V9 . (52) 
It can be readily shown that the maximal A along h = 29 is 

A max = -lmV2 + l) ~ 0.28 (53) 

7T 



and the average in the instability band is A ~ 0.22 . It is evident that in the large 9 
region the growth rate never diminishes, always with a sizable constant A. Thus the 
exponent of particle production rate grows roughly in the time interval of 1/A (~ 5) x 
oscillation period 

Particle production is halted by the back reaction. This problem is studied by 
solving a coupled system of the inflaton and the radiation energy densities in the 
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expanding universe p6| . 

+ 3Hn — N± M , (54) 
^ + AHp r = Nj t ( Ptp ) , (55) 

H 2 = yG(ft + fr), (56) 

where N is the number of boson species contributing to parametric resonance effects, 
roughly with the coupling g. 

The idea behind this evolution equation is that irrespective of interaction among 
created particles towards thermalization, the energy balance between the inflaton 
and radiation should hold at any instant of time because of energy conservation. If 
one separately checks that thermalization is realized, the radiation energy density p r 
can be used to estimate the temperature in equilibrium; 

We allow the possibility that there might be more particle species participating in 
equilibration than the number of created boson species N: for instance, in supersym- 
metric theories of N boson species there are roughly equal number of fermion species 
that may be produced from the bosons by secondary processes, giving N' ~ 2N. The 
system of evolution equation can be extended to include the mass variation by 

¥ = < 58 > 

It might be instructive to give derivation of this evolution equation in view of 
that the meaning of energy densities in the presence of interaction terms is not clear. 
We define the inflaton energy density including the mass variation; 

1 ^2 1 
-f + - 
2 s 2 

with the time dependent part of the mass given by g 2 Y^<p (y? 2 ) • From the £ equation 
of motion, 

f)V 

t+3HZ + ml{t)Z = -(-£), (60) 

where (A) = tr (Ap) (p = density matrix) and Vy = \ gm^ J2 V f 2 , one derives for 
evolution of the £ energy density; 

P, + 3^ = -(^>4^X> 2 >. (61) 



P 5 = ^ 2 +o"#K 2 ' (59) 
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Note that we included the quartic interaction term | g 2 £ 2 J2 V (v? 2 ) i n the time variant 
mass term m 2 {t). The left hand side becomes p$ + 3Hp^ after time average over the 
oscillation period, since 

e=^f(tye = P?, (62) 

under the assumption of slow £ variation. 

On the other hand, it can be shown using the time evolution equation for the 
density matrix, 

p=-i[H,p] = -i['£H ip + V Y + V A ,p], (63) 



that 



where 



tf, = ^ 2 + ^(V<p) 2 (65) 
is the free Hamiltonian density of created boson 99 and V± = \ g 2 ^ 2 J2 V <P 2 ■ The 



short time average over one oscillation period of (Vy) vanishes, since this quantity 
has a linear dependence on £. Thus one has 



leading to 

p f + 3^ = - ^ E W . ( 67 ) 

p r + 4tf p r = |E . (68) 

The second equation simply follows from energy conservation of the combined system 
of £ and radiation. 

It is now appropriate to discuss the thermal history after inflation, namely after 
the inflaton commences oscillation at t ~ § — . The result should be compared with 
the naive estimate of the reheat temperature, T B ~ 0.1 v /7^m p i ~ 1CT 2 g m v \ . 
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Particle production is significant only when the exponent in the rate is appre- 
ciable, but once it becomes of order unity, the production rate becomes accelerated. 
Indeed, when the exponent becomes of order 100, a catastrophic particle production 
occurs and back reaction immediately stops particle production. The time when this 
happens can be estimated by 

t = t d aO[100]/(Xmt). (69) 

This is an abrupt change, as numerically checked by solving the time evolution 
equation [3^]. Prior to this time the inflaton density varies according to the £— matter 
dominance, p^ = |m|£ 2 oc 1/t 2 . Combined with the time of catastrophic decay td 
above, the abrupt change occurs at 

JUo[KT 2 ]|a. (70) 

A typical energy of produced particles at this time, prior to thermalization, is of 
order 



E = O[ s Jgm^] = O[0.3]^g\m $ m pl , (71) 

with £ given as above. Using this energy, one estimates the rate of two body reactions 
among created particles, and the corresponding ratio of this to the Hubble rate, 

T_ ^ a s /E 2 ■ Nn v 
H sfppTpr '/m p i ' 

m pl / Pr s 1/2 



E -E 4 \l + Pr /pt 



Pr/Pt (?2) 



with a s a typical coupling involving ordinary particles. What usually happens |36| 
is that 

Pr^Pt, ^>1, (73) 

hence T ^> H . Thus two-body reactions frequently take place. Even multi-particle 
reactions occur and one may conclude that thermalization takes place immediately 
after the catastrophic particle production. Under this circumstance the reheat tem- 
perature right after the explosive particle production is 

T R ^O[10- 2 ] y/Xm^m pl . (74) 

With ~ 10 13 GeV , A ~ 0.2 , this reheat temperature is close to the GUT scale. 
It is possible that at least the Higgs X— boson of mass ~ 10 14 GeV can copiously be 
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produced. In any case T R 3> T B for a small coupling g. More precise computation 
has been done by numerically integrating time evolution equation . 



The entire thermal history may in general be fairly complicated. Even if the 
radiation dominance is realized at the early epoch simultaneous to the catastrophic 
production, there exists a residual inflaton field after the catastrophic stage: 



-4i 1 /^Ampi 2 



6 « O[10" 4 ] -C-^Y . (75) 
Six rri£ 

If this value is still large, ^> 1, there may be a second explosive decay. On the 
other hand, if this value is moderately small, but not small enough such that the 
naive perturbative analysis is no longer valid, then there may be gradual particle 
production continuously down to 9 <C 1. This region is difficult to analyze, but 
under study currently. 

The ultimate end point of the inflaton decay is the decay in the first band, 
namely the Born decay when the inflaton amplitude becomes very small. This is 
because only the first band appreciably contributes with very small 9, and the Born 
decay rate always satisfys 7g > H at very late times. Under the assumption that 
intermediate amplitude decay is not significant, one can estimate entropy creation at 
the Born decay. According to our numerical integration |3(J the entropy production 
by the Born decay is significant for the coupling range g < O[10~ 3 ], but negligible 
for g > O[10~ 2 ]. In all cases numerically checked ||36|| , the final temperature at 
the end of the complete inflaton decay is given approximately by the estimate due 
to the Born formula Tb- This does not mean that effects of parametric resonance 
can be forgotten. High energy processes that may occur in the temperature range, 
T R > T > T B , have to be reconsidered. As such, the baryogenesis and the gravitino 
production are of prime importance. 



5 Gravitino abundance 

The old formula of the gravitino abundance n^/2/s ~ 1CT 2 • Tg/m p i with s the 
entropy density is not valid since the thermal history right after the catastrophic 
particle production is complicated in some region of the coupling g, and Tb is not 
a good measure to characterize the thermal history. The truth is that no single 
temperature represents the state after inflation. In order to estimate the gravitino 
abundance with non-trivial thermal history after inflation, it is necessary to follow 
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time evolution of the gravitino number density n 3 / 2 [36] , 



d ^+3Hn 3/2 = (i:v)nl, (76) 



where n v is the thermal number density of one species of created bosons. The 
cross section (Ev) of gravitino production <pcp — > g 3 / 2 #3/2 has been computed, and 



roughly (£t>) ~ 250 /m^ unless the gravitino is lightest supersymmetric particle 
(LSP). Destructive term has been neglected in the evolution equation above, which 



is justified for the gravitino mass larger than 0[1] keV [37 



Both possibilities of stable and unstable gravitino remain viable |36[ . Let us only 
mention a possibility of the gravitino dominated universe at the present epoch [|37| . 
With the initial temperature Tr > 2 x 10 14 GeV imposed to give a favorable situation 
for GUT baryogenesis, there is a region of parameters for the closure density of 
gravitino dominated universe if 1713/2 = 0.1 — 10 GeV. The basic reason this becomes 
possible is that the initial large gravitino yield created right after the catastrophic 
particle production is much diluted via the late phase of Born decay. This is again 
a reflection of the large disparity of the two temperatures; T R ^> T B . Of course, it 
remains to demonstrate a sizable baryon to photon ratio. But things are not bad: 
there is an epoch immediately before the catastrophic particle production in which 
non-equilibrium environment necessary for baryon generation exists, and moreover 
the observed baryon to photon ratio is of order 10~ 10 , allowing some amount of 
dilution in later epochs. 
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